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Abstract 

Motivated by the interplay between structural and reduced form credit models, we pro- 
pose to model the firm value process as a time-changed Brownian motion that may include 
jumps and stochastic volatility effects, and to study the first passage problem for such pro- 
cesses. We are lead to consider modifying the standard first passage problem for stochastic 
processes to capitalize on this time change structure and find that the distribution functions 
of such "first passage times of the second kind" are efficiently computable in a wide range 
of useful examples. Thus this new notion of first passage can be used to define the time of 
default in generalized structural credit models. Formulas for defaultable bonds and credit 
default swaps are given that are both efficiently computable and lead to realistic spread 
curves. Finally, we show that by treating joint firm value processes as dependent time 
changes of independent Brownian motions, one can obtain multifirm credit models with 
rich and plausible dynamics and enjoying the possibility of efficient valuation of portfolio 
credit derivatives. 

Key words: Credit risk, structural credit model, time change. Levy process, first passage time, 
default probability, credit derivative. 
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1 Introduction 



The structural approach to credit modeling, beginning with the works of Merton [|29ll and Black 
and Cox [|71, treats debt and equity as contingent claims (analogous to barrier options) on the 
firm's asset value process. This unification of debt with equity is conceptually satisfying, but in 
practice the approach sometimes leads to inconsistencies with intuition and observation, such 
as the zero short- spread property (a consequence of the predictable nature of the default time), 
and, in Merton type models, time inconsistency. The Black-Cox framework, while time con- 
sistent, also leads to technical difficulties when pushed to provide realistic correlations between 
different firms' defaults and with other market observables. Both approaches tend to be too 
rigid to allow good fit to market data as well as effective hedging strategies. 

Reduced-form (or "intensity-based") modeling, introduced by Jarrow and TurnbuU |!2T1|, has 
been successful in providing remedies for some of these problematic aspects. It treats default 
as locally unpredictable, with an instantaneous hazard rate, but does away with the connection 
between default and the firm's asset value process. 

Subsequent developments have bridged the gap between reduced form and structural mod- 
els. For example, the model of Jarrow, Lando and TurnbuU [19] and its extensions [|23l [Tl[T7l 
posit a continuous time Markov chain to replace the firm value process as a determinant of 
credit quality, while retaining the concept of hazard rate in the form of dynamically varying 
Markov transition rates. The time of default is the first-hitting time of the default state, an ab- 
sorbing state of the Markov chain. The incomplete information approach of Jarrow and Protter 
||20|| views reduced form methods as arising from structural models when the market has less 
than perfect information about the firm value. So-called hybrid models [|27ll (see also [11]) seek 
to tighten the connection with structural models by allowing the hazard rate to depend on the 
firm's equity value (stock price), and allow the stock price to jump to zero at the time of default. 

Another useful way to generalize the structural framework is to consider other classes of 
stochastic processes for the firm value. This line of work can potentially parallel a great body 
of work that extends stock price return models by allowing jumps and stochastic volatility. 
While this can be quite successful in the Merton framework, the Black-Cox framework faces 
the technical difficulty associated with the first passage problem. Black-Cox style structural 
credit models using jump-diffusion processes to model the firm value have been studied in 
ll32l [3TII . Because of the difficulty in solving the first passage problem, their models are hard 
to compute, but these studies do demonstrate the viability of the approach in curing some of 
the deficiencies of the classic Black-Cox model by adding flexibility and the possibility of 
unpredictable defaults. As for exact formulas for first passage, results based on fluctuation 
theory [6] and Wiener-Hopf factorization f5l are known. Kou and Wang [22] solve the first 
passage problem for a specific class of jump-diffusion process, and Chen and Kou [12] use 
those results to extend the Black-Cox firm value model and the Leland-Toft model [|25l for the 
optimal capital structure of the firm. 

The purpose of the present paper is to explore time-changed Brownian motions (TCBM) for 
their potential to be used for consistent modeling of a firm's asset value process and the firm's 
time of default. We aim to retain flexibility (to be able to match a wide range of possible credit 
spread curves), computational tractability (to permit efficient option valuation), and logical con- 
sistency with the paper of Black and Cox (by treating default as a first passage time for the firm 
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value to hit a default threshold). 

Many authors have used TCBMs as models of log stock returns (a notable review paper is 
ifTSll ). and their great flexibility is by now well known. When the time change is an independent 
Levy process (Levy subordinator), one obtains well known models such as the variance gamma 
(VG) model and the normal inverse Gaussian (NIG) model. Barndorff-Nielsen and Shephard 
ill have introduced time change models where the time change is an integrated mean-reverting 
jump process, while important stochastic volatility models such as Heston's model [il6j arise 
from time changes that are integrated mean-reverting diffusions. The paper of [fTOll surveys 
25 different realizations of time changed Levy processes and shows how they perform in a 
calibration exercise to observed option pricing data. 

This paper focusses on the remaining difficulty in successfully implementing TCBMs to 
model the firm value in a Black-Cox framework, namely the efficient computation of first pas- 
sage probabilities. Because of the difficulties that arise in computing the associated first passage 
distribution, and in analogy to the time-changed Markov chain models where the default state is 
an absorbing state, we are lead to propose a specific variation of first passage time applicable to 
time-changed Brownian motions, but not to general jump diffusions. This variation, which we 
call the first passage time of the second kind, is designed to be decomposable by iterated condi- 
tional expectation, and thus can be computed efficiently in cases of interest. This concept is not 
new, having been used for example by Moosbrucker [|30ll and Baxter p| in their computations 
of basket credit derivatives, but to our knowledge its modeling implications have not yet been 
fully explored. 

Our purpose here is threefold. First we explore the mathematical structure of first passage 
times for time-changed Brownian motion, and provide a set of natural solvable examples that 
can be used in finance. By comparison of these examples with a range of existing stock price 
models, we thereby demonstrate the broad applicability of our framework to equity and credit 
modeling. Our second aim is to focus on structural models of credit where the firm value process 
is a general time-changed Brownian motion and the time of default is a first passage time of the 
second kind. We prove pricing formulas for defaultable zero coupon bonds and credit default 
swaps, with and without stochastic recovery. This discussion demonstrates that time-changed 
Brownian motion can be the basis of single firm credit models consistent with the principles of 
no arbitrage, and with tractable valuation formulas for important derivative securities. Finally, 
we demonstrate how the single firm model can be extended to the joint default dynamics of 
many firms. Under a restrictive assumption on the correlation structure, analogous to the one- 
factor default correlation structure in copula models, we demonstrate the efficiency of valuation 
formulas for portfolio credit derivatives. 

To avoid obscuring our most important results by focussing on a too-specific application, 
we ask the reader's indulgence to postpone statistical work on the modeling framework to sub- 
sequent papers. While we are hopeful that statistical verification of the modeling assumptions 
on asset price datasets will ultimately show the viability of our framework, such a verification 
must proceed one application at a time, and would take us too far in the present paper. 

In outline, the paper proceeds as follows. Section 2 introduces the probabilistic setting 
and the definition and basic properties of TCBMs. The first passage problem for TCBMs is 
addressed in Section 3. Since the standard first passage problem for TCBMs exhibits no sim- 
plification over first passage problem for general jump-diffusions, we introduce an alternative 
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notion, called the first passage time of the second kind, that capitalizes on the time change struc- 
ture. It is this notion that is used in all subsequent developments. Section 4 introduces the main 
categories of time changes, namely the Levy subordinators and the integrated mean-reverting 
jump-diffusions. These two families are in a sense complementary, and together provide a rich 
and tractable family of TCBMs that can be used as building blocks for more general models. 
Section 5 introduces the simplest structural credit models based on TCBMs, and runs through 
the valuation of some basic credit derivatives. Section 6 provides a brief numerical exploration 
of the single firm model. The multifirm extension is addressed in Section 7. We find that com- 
putational tractability strongly suggests that while the time change processes for different firms 
may (indeed should) be correlated, the underlying Brownian motions must be taken independent 
firm by firm. 

2 Time-changed Brownian motion 

Let {n, JF, F, P) be a filtered probability space that supports a Brownian motion W and an 
independent strictly increasing cadlag process G with Gq = 0, called the time-change. P may 
be thought of as either the physical or risk-neutral measure. Let Xt = x + aWt + (3a'^t be the 
Brownian motion starting at x having constant drift Pa'^ and volatility a > 0. 

Definition 1. The time-changed Brownian motion (TCBM) generated by X and G is defined to 
be the process 

Lt = XG,, t>0. (2.1) 
To eliminate a redundant parameter, we normalize the speed of the time change: 

lim T-^E[Gt] = 1. 

T — ^■OO 

In addition to the "natural" filtration of the TCBM, J^t = a{Xs,Gu ■ s < Gt,u < t}, 
one can also consider the subfiltration Qt = a{Gs ■ s < t} and the Brownian filtration Wt = 

a{Ws : s<t}. 

Remarks 2. 1 . Since the above definition requires X and G to be independent, it is some- 
what more restrictive than notions of TCBM studied by others. However, as has been 
amply demonstrated by finance researchers such as [fTSlfTOl . this family of processes of- 
fers a promising degree of versatility and tractibility when used to model equities. We 
will see in what follows how this family can be useful in credit risk modeling. 

2. It is a point of philosophy to think of the time change as a reflection of the impact of 
the market on individual firms, in the same way that in reduced form models, the default 
intensity A lies in the market filtration. We think of the processes X as firm specific 
generalizations of the default indicator function, while the time change generalizes the 
default intensity. 
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The characteristic functions $ for any Q < s <t and u e V,V a domain in C are defined 
to be 

^[gm(it-i.)|jr^]. (2.2) 

and real variable versions of these functions, in the form of the Laplace exponents ip{u) ~ 
— log$(iw), are also of interest. All are to be understood as processes in the variable s. A 
simple calculation gives an essential formula 

= E[^l{u,Gt)\gs] = ^^{a'{Pu + zu'/2),t). (2.3) 

This formula reflects a kind of "doubly stochastic" or "reduced form" property of TCBMs, 
already familiar in credit risk modeling: general JF expectations involving default can be "re- 
duced" to expectations in the "market filtration" Qg. 

Useful "solvable models" arise when $f and hence are explicit deterministic functions 
of an underlying set of Markovian variables. Explicit characteristic functions not only lead to 
formulas for moments m^''^ = ElL'l] and cumulants c*^'^^ for k = 1,2, . . . but are also useful for 
extracting more detailed properties of the process L. An important algebraic aspect of TCBM 
is their natural composition rules. If G,H are two independent time changes then G + H, 
G X H and G o H are also TCBMs, and one has results such as $^+^ ^ x and 



3 First passage distributions 

The first passage time for a semimartingale Lf starting at a point Lq = x > to hit zero is the 
stopping time 

t(^) = inf{t|Lt < 0}. (3.1) 

The first passage problem for general semimartingales is difficult to deal with, and few appli- 
cations have been realized. When applied to the nice subclass of TCBMs, it turns out that first 
passage times do not respect the reduced form property, and even in this setting they remain 
difficult to compute. On the contrary, when Lf is a TCBM process, the following different def- 
inition of first passage time, which we call the first passage time of the second kind, is both 
natural and compatible with the "reduced form" property. Consequently we shall see that it can 
be easily implemented in credit modeling. 

Definition 3. The first passage time of the second kind of the TCBM Lt — Xq^ is the stopping 
time 

^(2) = mi{t\Gt > t*}, (3.2) 

where t* = mi{t\Xt < 0}. 
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Remarks 4. 1. As shown in [18], t^^^ can be viewed as an approximation of the usual first 
passage time t^^^ with t^^^ > t^'^\ When G is a continuous process, the two definitions 
coincide. 

2. When the time change is a pure jump process with unpredictable jumps, both stopping 
times are totally inaccessible. In general, they can be written as the minimum of a pre- 
dictable stopping time and a totally inaccessible stopping time. 

An essential first step in constructing tools for studying t^"^^ is to collect "structure" functions 
associated with first passage t* for drifting Brownian motion Xt = x + aWt + Pa'^t. The 
following formulas are well known (see for example [91): 

1 . The cumulative distribution function for the first passage time of drifting Brownian mo- 
tion is 

Pit, X, a, P) := i?.[l|,.<,|] = N (^^^^) + e-'^^N (^^^^) • (3-3) 

2. For any u > — /?^a^/2, the Laplace exponent of t* is 

il:{u,x,a,P) := -logE^[e-"**l{j*<oo}] 



-log 



X 



3. The joint distribution function Ex[l{f>t}'^{Xt>e}] is given for £ > by 

These elegant formulas of Brownian motion are needed in the theory of the second kind of 
passage problem, for which the structure functions of t^"^^ are computable via an intermediate 
conditioning. Thus, for example, its cumulative distribution function (CDF) is given by: 

P^^\t,x) := E,[l|,(2)<,}] =E[E,[l|i*<G,}|6?oo]] (3.6) 
P{y,x,a,f3)pt{y)dy 



where pt is the density of Gt and the function P is given by p.3| ). 

While this formula can sometimes be used directly, in the many cases of interest where the 
Laplace exponent ^(m, t) of the time change Gt is given in closed form, P*^^^ and other structure 
functions can be given a more useful Fourier representation. 

Proposition 1. For any x > let Lt = Xct ,Xt = x + aWt + /Scr^t be a TCBM where Gt has 
Laplace exponent ■ip{u,t) := — logi?[e~"'^*]. Then 
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1. For any t > and e G M the function i?^[l|j^((2)}(5(Lt — £)] is given by 



4^>o}- 



27r 



(3.7) 



while Ex[^!^y^{2)^5{Lt — ^)] is given by 



27r 



(3.8) 



M+ie 



2. For any t > ^/ze CDF P*^^) (t, x) ^/ven 

e i/3>0 + -'-/3<0 — 



e '^^ zsmizx) 



TT 



z2 + /52 



(3.9) 



while the characteristic function Ex[li^^^^{2)-^e is given for any k in the upper 

half plane by 



-I3x 



271 



k — z k + z 



(3.10) 



Remark 5. The formulas in this proposition are all explicit Fourier integrals in the variable x 



involving the Laplace exponent. This is a key advantage over a formula like p.6[ ) in at least two 
respects. Firstly, the Laplace exponent is in many cases given explicitly while the density p is 
not. Secondly, compared to a generic numerical integration, the fast Fourier transform (FFT) 
algorithm is computationally efficient and comes with powerful error estimates as described in 

m. 



Proof of Proposition [T]: To show p.7[ ), we note that 



where the inner expectation can be evaluated by differentiating ( |3.5| ) and using a standard Gaus- 
sian integral that holds for any e and t > 0: 



This then leads to 



27r 



-izx—z^a'^t/2^^ 



(3.11) 



E 



Ml-x)-l3^a^Gt/2 



2n 



-{£>o}- 



-{£>o}- 



Mi-x) 

2ti 

,j3{e-x) 

271 



^ize ^^-izx _ ^izxi^ E[e-"'("'+^')^*/2]^2 



izl 



where in the second step we have used the Fubini Theorem to interchange the integral and 
expectation. 

To prove p.8[) we use similar logic to note that 



E,MLt-^)] = E 



1 

2^ 



To prove p.9[ ) we integrate ( |3.7| ) over ^ G To interchange the order of the ^ and z 
integrations we need to take e > /5 if /5 > 0, but may take e = if < 0. Then we find 



-Px 



X) 



2-K 



-I3x 



27r 



1 



'y^-izx _ ^izx^ g-V'(<T2(^2+/32)/2,t)^^ 



i{z — i(3) 



The integrand has a residue of (27ri)~^e~^^'* at z = ijS, and when /3 > we need to apply the 
Residue theorem to deform the ^-contour to the real axis. For all (3, the resulting integrals over 
M can be manipulated into the required real form. The proof of p.lO| ) is straightforward under 
the condition on A;. □ 



4 Time change models 

TCBMs have been well studied as models of log-stock prices. We outline a general approach 
to building TCBMs, and provide a number of distinct types of time change that can be used as 
"building blocks". 

4.1 Levy subordinated Brownian motions 

These TCBMs arise by taking G to be a Levy time change, that is, a Levy subordinator. Such 
TCBMs are examples of Levy processes, the general class of continuous time stochastic pro- 
cesses with stationary and independent increments. Much of the analysis connected with a Levy 
process Lt is based on its characteristic triple (6, c, p)h, in terms of which its Laplace exponent 
takes the form 



ilj^{u,t) := -logE[e-"^*] = t 



bu - c'u^jl + / [e""^ - 1 - uyh{y)]p{y)dy 



(4.1) 



Here p is a measure on M \ 0. For ease of exposition in what follows, we set the truncation 
function h{y) to zero, which is permissible by adopting the restrictive condition that |x| A 1 
should be p-integrable. Our main results extend to the general case where |xp Al is p-integrable. 
See [13J for general discussions of Levy processes. 

The following result is Theorem 4.3 in [13|, and identifies the type of process that can be 
expressed as a Levy-subordinated Brownian motion (LSBM) Lt := Xg^'- 
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Theorem 2. Supposing Lq = x, the following are equivalent statements: 

1. L is a Levy process with characteristic triple (6, c, p)o where c > 0. The density p is 
nowhere zero on M and for some 13, e^v^p(-v^) = e-^v^p(y^) and is a completely 
monotone function. Furthermore h = (3c. 

2. Lf := for drifting Brownian motion Xt = x + Wt + (3t and G a Levy subordinator 
with characteristic triple {b, 0, z/)o, b > 0, and v a measure on (0, oo). 

Here are some examples time changes G for Lt = Xc^iXt = x + aWt + (3aH that have 
been used in models of logarithmic stock returns: 

1. The exponential model with parameters (a, 6, c) arises by taking G to be the increasing 
process with drift 6 > and jump measure v{z) = ace~"'^,c,a > on (0,oo). The 
Laplace exponent of Gt is 

^{u,t) := - \ogE[e-''°'] = t[bu + uc/{a + u)]. 

and the normalization condition is b + c/a = 1. The resulting time-changed process 
Lt := Xct has triple {/3b, b, p)o with 

where (?/)+ = max(0, y), (?/)"" = (— y)^. This forms a four dimensional subclass of the 
six-dimensional family of exponential jump diffusions applied to finance in [|22ll . 

2. The VG model ll26l arises by taking G to be a gamma process with drift defined by the 
characteristic triple (6, 0, u)q with 6 > (often b is taken to be 0) and jump measure 
^{z) = ce~"^/z, c, a > on (0, oo). The Laplace exponent of Gt is 

^{u,t) := -logE[e""^*] = t[bu + c\og{l + u/a)]. 

and the normalization condition is6 + c/a=L The resulting time-changed process has 
triple {/3b, b, p)o with 

\y\ 

3. The normal inverse Gaussian model (NIG) with parameters /3, 7 O arises when Gt is the 
first passage time for a Brownian motion with drift /? > to exceed the level ^t. Then 

^{u,t) =^t{/3+ ^J/3^ + 2u) 

and the normalization condition is 7//? = 1. The resulting time-changed process has 
Laplace exponent 

il^^{u,t) =xu + t^[/3+ J (3"^ + 2(3u + V?]. 
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4.2 Affine TCBMs 

For our second important class of time changes, Gt = Jq Xsds has differentiable paths, and the 
corresponding TCBMs are diffusions (processes with continuous paths) which exhibit "stochas- 
tic volatility". We focus here on a class we call ATCBMs ("affine" TCBMs), for which A is 
taken in the class of positive mean-reverting CIR-jump processes introduced by [fT4l . We men- 
tion here two distinct examples: 



dx[^^ = {a-bX^^^)dt+ \/2cXi^^dW}^\a,b,c>0, 

dXf^ = -hX^^^dt + dJt. (4.2) 

Here J is taken identical to the exponential Levy subordinator with parameters (a, 0, c) defined 
in example 1 of the previous subsection. 

The essential computations for Laplace exponents 

V^(*)(u,t;A) := -logE[e-"^*''|Aj^ = A],i = 1,2 

of such affine time changes are described in many papers. The following formulas are proved 
in the appendix of f[T7l : 

Proposition 3. The characteristic functions ip'^'^\i = 1, 2, both have the affine form 

^^'\u,t-X) = A»(M,t) + A5«(M,t). (4.3) 
The functions A^*^ and B^'^^ are explicit: 
1. 

A«(M,t) = -K2 + (l + (e^* - l))^'/€2, 

5(1) (u, t) = -aKit + I log {\ + ^/ti (e^* - 1) 
with constants Ki, K2 cmd 7 given by 



(4.4) 




(4.5) 



A(2)(M,t) = I (^1-e-^*) , 

t) = ct-^ log ( ^^^+^)f~^ 

^ ' ab+u \ 0,0 



(4.6) 



The ATCBM model with time change X^J^^ds is equivalent to the Heston stochastic volatil- 
ity model for stock returns fil6il . with zero correlation (hence zero leverage effect). Stock price 
models with time change /* Xf'^'ds, and extensions thereof, were introduced in [|2l. 
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4.3 More general TCBMs 



Two different ways of combining time changes have been studied that preserve the desirable 
property that the resulting Laplace exponent is explicitly known. The first is to add time 
changes. For example, a model that includes time jumps and stochastic volatility with both 
a diffusive and jump component arises if we take 

G,= [\xW + Xf))ds + G['^ 
Jo 

where Gf'^ is a Levy subordinator. 

The paper [[T0| explores 25 models of the form Lnt where L is a Levy process and Ht = 
Jq Xsds is an independent time change similar to those mentioned in this section. In most cases 
they discussed, L is itself of the form Xct, and so is a TCBM with time change G o H. 
Amongst this class, the paper finds several models of the stock price that capture very well the 
implied equity volatility surface. 

From the large literature on such TCBMs as models for equities, it is clear that the class of 
TCBMs with explicit Laplace exponent is rich enough to describe a wide range of asset classes 
in finance. We shall now see how such processes can be used to build structural models of credit 
risk. 



5 Structural credit models 

The structural credit modeling paradigm of Black and Cox [|7l assumes that default of a firm is 
triggered as the debt holders exercise a "safety covenant" when the value of the firm falls to a 
specified level. It makes sense therefore to assume that the time of default is the time of first 
passage of the firm value process Vt below a specified lower threshold function K{t). 

In this section, we outline how a Black-Cox credit framework can be built under the assump- 
tion that the log-leverage ratio Lt = \og{Vt/ K(t)) is a TCBM. To demonstrate the flexibility 
of the approach, we make use of all the building blocks introduced so far, leading to a large 
number of parameters. A more realistic implementation would likely begin with a much more 
restricted specification. In analogy to the multifactor reduced form modeling framework of [|T4l| 
we choose the approach in which independent time changes are combined together by addition 
rather than composition. The alternative route via composition is deserving of separate study. 

Assumptions 1. 1. There is a vector Zt = [r^, X[^\ X^f''] of independent processes with A^*^ 



chosen as in subsection 4.2 r is a CIR process with Laplace exponent ilj'^{u, t) given in 
the form (|44j). 



The log-leverage process Lt = \og{Vt/K{t)) = Xa^Xt = x + aWt + (3aH is a TCBM 
where the time change is a convex combination of the building blocks of the previous 
section: 

Gt = a^Gl^^ + aaCf ^ + agCf ^ (5.1) 

with < tti, a2, as < 1 = «! + 02 + as- Here G[^^ = /J X^s^ds, i = 1,2 are defined as in 
Section 



4.2 



with Laplace exponents tp^'^\u, t; A*^*^) while G^^^ is a Levy subordinator with 



Laplace exponent ijj^^\u, t). We also assume that P < 0, a > 0. 
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3. The time of default is t'^^\ the first passage time of the second kind. 

4. The spot interest rate is n = rt + rriiXl^^ + m2\f'^ for non-negative coefficients mi, 1712. 

5. A constant recovery fraction R < 1 under the recovery of treasury mechanism is paid on 
defaultable bonds at the time of default. (This is for simplicity: [fTTll shows that we can 
allow Rt to be a general affine process. We can also compute under the recovery of par 
assumption with somewhat more complex integrations.) 

As previously mentioned, Black-Cox models using jump diffusions are usually based on the 
standard first passage time, leading to technical difficulties that can be solved in only a restricted 
class of processes. Our innovation is to consider instead the second kind of first passage time, 
and thereby capitalize on a type of "reduced form" for computations that applies to general 
TCBMs. 



5.1 Bond pricing 

The following proposition gives formulas for default probabilities and default-free and default- 
able zero coupon bond prices. 

Proposition 4. Let the initial credit state of the firm be specified by initial values Lq = x and 
= [fo, \f\ Recall that /? < 0. 

1. The probability that default occurs before t > is given by 



P[t(2) <t] = l 



-fix fOO 



^exp 



z2 + /52 



-J]^«(a.a2(^2 + /32)/2,t) 



dz. (5.2) 



2. The time price Pq (T) of the default-free zero coupon bond with maturity T is 



Po(T) = exp 



j=l,2 



(5.3) 



3. The time price of the defaultable zero coupon bond with maturity T, under constant 
fractional recovery of treasury, is Pq^{T) = (1 — R)Pq{T) + RPq{T), where Pq(T) 
denotes the price of the zero recovery defaultable zero coupon bond, given by 



z sm[zxi 



X exp 



=1,2 



(5.4) 



dz. 
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Proof. We prove only the formula for Pq{T): the other formulas are similar, but easier. The 
risk-neutral pricing formula gives 



P°(T) 



E 
E 



'e-!^^^'^^E[l{G,^,,}\g] 



-(3x 



Z ^Yil\ZX] 



E 



dz 



The third equality comes by using p.9[ ) in the case when Gt is deterministic, and using the 
Fubini Theorem to interchange the order of integration. The final form for ( |5.4| ) now follows by 
decomposing r^ds and Gr into their independent components and performing the resulting 
one dimensional expectations. □ 



The above pricing formulas are explicit functions of the initial values Lq = x, Zq; as time 
develops, prices are deterministic functions of the processes Lt and Z^. We adopt the point of 
view that Z contains information about the drivers of general credit markets, while L reflects 
firm specific information. 



5.2 Credit default swaps 

We next consider an idealized CDS with unit notional and maturity T > 0. The premium leg, 
paid continuously at a constant rate S until t*^^) /\ has the time price 



S V{T) = SE 



S I P^{t)dt 
'o 



(5.5) 



where the second equality comes by comparison to ( |5.4| ). The default leg pays the fractional 
loss of treasury value (1 — R)Po{t^'^^) at the time of default if t'^^) < t and has the time value 
W{T) = i^(Po(T) - Po°(^))- The CDS spread S{T) is defined to be the value of S that 
makes SV = W. 



6 Numerical results 

The structural credit modeling framework of the previous sections is designed with flexibility 
and computability in mind. Rather than embark here on a lengthy statistical investigation of 
promising specifications and their calibration to market data, instead, in this section we strip 
out the complexity, and simply exhibit a set of parametrizations of the VG TCBM model that 
generate plausible credit spread curves, thereby demonstrating the computational efficiency. 
We consider the credit framework for a pure geometric Brownian motion model (Model A) 



and three parametrizations of the VG model of ^: 4.1 (Models B,C,D), with parameters shown 
in Table 1. All four models are specified so that Lt has Lq = 1.5, fixed annualized variance 
+ /3V^(2/a + 1) = 0.09 and mean log rate of return ~a^/2 (i.e. /3 = -0.5). 
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In Figure 1 , we compare the thirty year zero recovery yield spread and defauk probability 
density for these models. We observe that the yield spreads equalize as maturity increases, but 
show the completely different short time behaviour expected from the presence of jumps. Figure 
2 shows the thirty year zero recovery yield spreads in Model B for four firms which differ in 
their initial distance-to-default values Lq = 0.3, 0.6, 1.0, 2.0. We see here that firms with small 
Lq (high default risk) can have decreasing spread curves, while the reverse is true if Lq is large. 





Model A 


Model B 


Model C 


Model D 




1.5 


1.5 


1.5 


1.5 


a 


1 


1 


10 


100 


b 


1 











c 





1 


10 


100 


(3 


-0.5 


-0.5 


-0.5 


-0.5 




0.09 


0.0846 


0.0877 


0.0880 



Table 1 : Parameter values for the VG TCBM model. 



7 A one factor multifirm structural model 

A difficulty in structural credit risk modeling is finding a framework that extends naturally and 
efficiently to a large number of firms, while allowing for a rich default dependence structure. 
The present setup of time-changed Brownian motions is such a framework. Consider M firms, 
where for each m = 1, 2, . . . , M, the mth firm is governed by its firm value process V^'^, default 
trigger threshold K™'{t) and log-leverage ratio process 

= logV^'^/K'^it) = X^^, 
Xr = x^ + a^Wr + l3^alt. (7.1) 

For the mth firm we take parameters Xm, /?m, o"m > 0, and a firm dependent time change G™. 
Assumptions 2. The joint dynamics of multifirm defaults is determined by the first passage 

(2) 

times tm of the log-leverage ratio processes L™. The time change processes G™ are given 
jointly in the "one-factor" form 

GT = a^Gt + (1 - a^)H^ (7.2) 



with am € [0, 1] and time changes G, if™ having the form given by ( |5.1[ ). Finally, we assume 
that G, if^, . . . , if^^, X^, . . . , X^^ are mutually independent processes. 

In this one-factor time change model, the maximal correlation structure is obtained by set- 
ting each am = 1- However, since the underlying Brownian motions X™ are independent, 
maximal correlation does not mean the defaults are fully correlated. 

The one factor model can be interpreted as a generalized Bernoulli mixing model, in the 
sense of [8] and [|28.I , where the mixing random variable is Gt- That is, the default states 
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of all firms at time t are conditionally independent Bernoulli random variables, conditioned 
on Qt '■= crlGs : s < t}. If we define the conditional survival probability F"\xm,Gt) := 
Exm\^{Gt<t"'*}\Gt\^ then the following formula extends p.9| ) and is proved exactly the same 
way: 

F^{x,y)= {13) 
' "^^"""^ exp[-a„,a^(z^ + Pl)y/2 - ^^1(1 - a^)al{z^ + /3^)/2, y)]dz. 



J-oo + 13, 



Now, for any subset a C {1,2,..., M }, the unconditional probability that the firms in 
default at time t are precisely the firms in a is given by 

POD 

I J] (1 - y)) n F'^ix"', y)Pt{y)dy, (7.4) 

where pt is the distribution function of Gt- 

There are by now well-known techniques that under the assumption of conditionally depen- 
dent defaults, reduce the computation of credit portfolio loss distributions and CDO tranches to 
intensive computation of the conditional survival probabilities F™-{x, y). 



8 Conclusions 

We have studied the first passage problem for a class of semimartingales that are important for 
financial modeling, namely Brownian motions time-changed by an independent time change 
process. It was seen that the first passage time of the second kind presents some key advantages 
over the standard definition of first passage time, particularly computational tractability and the 
possibility of extension to multi-dimensional processes. 

Based on these good properties, we defined a pure first passage structural model of de- 
fault, and obtained computable formulas for the basic credit instruments, namely bonds and 
CDSs. The resultant formulas resolve a fundamental deficiency of the classic Black-Cox for- 
mula, namely the zero short spread property, and provide needed flexibility to match details of 
yield spreads. 

Finally, we outlined an extension to many firms in which dependence stems from systemic 
components to the time change, while the underlying Brownian motions are independent and 
firm specific. The resulting multifirm framework has a conditional independence structure that 
enables semianalytic computations of large scale basket portfolio products such as CDOs. The 
present paper focussed entirely on the mathematical properties of this modeling approach, and 
leaves interesting implementation questions such as calibration and applications in portfolio 
credit VaR and CDO pricing as subjects for future work. 

Acknowledgments: I am grateful to Alexey Kuznetsov and Zhuowei Zhuo for helpful discus- 
sions on the general theory of time changed Brownian motions. 
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Figure 1: Thirty year yield spread and default PDF curves for geometric Brownian motion 
model and three versions of the VG TCBM credit risk model. 




Figure 2: Thirty year yield spread for Model B with four different values Lq = 0.3, 0.6, 1.0, 2.0. 
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